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It is found that the problem of existence of bound entangled states with non-positive partial transpose (NPPT) 
has the intriguing relation to the Hilbert's 17th problem. More precisely, we compute the expectation value of 
the partially transposed Werner states by Schmidt rank-2 vectors for NPPT and 1-copy undistillable region. It 
is the positive polynomial but shown not to be expressed as the sum of squares of polynomials. A remedy for 
such pathological behavior as well as a similar but different approach to the problem is also mentioned. 
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O ' I. INTRODUCTION 

^_^ quantum information processing [1-3]. These fascinating tasks can be caiTied out with the highest performance by exploiting 

^ maximum entangled pure states. The pure entangled states are, however, very fragile and easily degraded to be mixed due to the 

'""I unavoidable interaction with the environment, whereas, in a quantum communication setting, two or more distant parties are in 

^NJ principle not allowed to perform global operations on the shared entangled states; their ability is quite limited even though they 

^-H , can communicate each other As a remedy for such an unwanted decoherence effect of entanglement within the limits of distant 

party setting, the concept and protocols of distillation of entanglement, the procedure to extract pure maximally entangled states 

(— ! . out of several copies of mixed states by using only local operations and classical communication (LOCC), have been introduced 

Q_i' |3t31- At presentjthe distillability — the possibility of distillation — is in turn recognized as one of the most fundamental traits 

^ ' of entanglement lIsl-ITlll. 

^ The distillability II12I1 of bipartite states is provided by the following theorem. A bipartite state p acting on a composite Hilbert 

C^ space "Ha ® "Hb is distillable if and only if there exist a positive integer A^ € N = {1,2,...) and a Schmidt rank-2 state vector 

^ \iJ/[^^ ) in •7/«^ ® -Hf such that CI 






(^ ■ Here, Tb denotes the partial transpose with respect to system B and the Schmidt rank-2 state vector is of the form, 

T-H ' ,=1,2 

J^ with Ci e R and |e,)^ (1/;)^) (/ = 1 , 2) being orthonormal vectors in 'K®^ CHf'^). 

f^ ' If Eq. ([T]i is satisfied for some A^, p is called A-copy distillable. Conversely, if it does not hold for up to some finite N eN, p 

T— I is called A-copy undistillable. While A-copy distillable states are distillable by definition, the A-copy undistillability does not 

►^ necessarily mean undistillability. To show the undistillability, we need to prove the A-copy undistillability for all A € N. 

• I— I . Not all entangled states are distillable. Undistillable and yet entangled states are called bound entangled (BE) states 111811 . One 

/^ ■ notable example of BE states is a state p whose partial transpose with respect to system A(B) remains positive semi-definite; 

?H p^'W*) > 1 18]. Such a state is called a positive partial transpose (PPT) state. Otherwise, it is called a non-PPT (NPPT) state. 

■ - - ■ It is clear from the above theorem that all PPT states are undistillable and that the distillable states are inevitably NPPT. All 

distillable states known to date belong to the NPPT class. So the natural question is: Is the NPPT condition also sufficient for 

distillability? or is there an NPPT BE entangled state? Since the NPPT condition is proven to be equivalent to distillability for 

two-qubit systems 1 19] and for states on "Ha ® "Mb with dim'TY^ = 2 and dim 'Kb > 3 [20], an NPPT BE entangled state, if it 

exists, must live in higher dimensions [21]. This is a long-standing open problem and remains unanswered since the discovery of 

bound entanglement in 1998 ll9l 4llll . Although considerable efforts have been devoted in solving this problem so far ll20ll22l - l33ll . 

the problem has been defying solution stubbornly. 
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Two early independent attempts ll20ll22|] (see also fls]) strongly suggest that NPPT Werner states lying at a narrow region 
adjacent to the PPT-NPPT borderline are undistillable; they have identified A^-copy undistillable NPPT Werner states but the A^- 
copy undistillable region is getting narrower and narrower and vanishes eventually as A^ goes to infinity. Although their findings 
may support the existence of NPPT BE Werner states, it is still inconclusive. Guessing from the finite A^ results may be quite 
misleading because there are states which are not A^-co py d istillable but (A^ + l)-copy distillable 1341 . The same is very much 
true for small A^ numerical study 12911 . A Werner state IBSll in question acting on ^Ha ® "Hb with dim'TY^ - dim "Kb = di> 3) 
takes the form, 

Pw(a) = 2_^ ^^ ~ °'^^' 
where F - 2f /=i \U) (i'l denotes the swap or flip operator and -1 < a < 1. The partial transposition of pw(a) reads 

where P+ = I^F+X^F+I with |*I'+> = d'^'^ Xti l">- For -1 ^ « ^ l/dpw{a) is PPT and also separable, while for l/d < a < I 
it is NPPT. The NPPT region is further divided into two; for 1/2 < a < 1 pw((^) is 1-copy distillable but for l/d < a < 1/2 it 
is 1-copy undistillable, i.e., there does not exist a Schmidt rank-2 vector |(/'2 ) ^ "^a ^'Hb such that uj/^ \pw{oY'' \4'2l < in 
this range of parameter a. The conjecture is that for l/d < a < 1/2 pw{o:) is also A^-copy undistiUable for all A^ e N, i.e., it is 
undistillable in an exact sense. 

It is pointed out that any NPPT state is distillable if and only if any entangled Werner state is distillable ll36ll . The immediate 
consequence is as follows: if NPPT BE Werner states do not exist, then the NPPT condition is necessary and suflicient for 
distillability. Thus, the solution would finalize the classification of (bipartite) quantum states from the viewpoint of entanglement 
whether the answer to the conjecture is positive or negative. Furthermore, the existence of NPPT BE states has the following 
important implications |37]. Namely, if they exist, then the distillable entanglement |5] — one of the entanglement measures — 
is not convex, and more surprisingly, the compelling yet counterintuitive phenomenon called superactivation [38] is expected, 
that is, the tensor product of two BE states, one some PPT BE state and another the conjectured NPPT BE Werner state, turns 
out to be distillable. It is furthermore shown in |41] that any NPPT state can be boosted to be 1-copy distillable by adding (tensor 
multiplying, to be exact) a PPT BE state with infinitesimal amount of entanglement. This implies that, if the conjectured NPPT 
BE exists, then the superactivation is always triggered by an infinitesimally small amount of PPT BE. 

In this paper, we point out that the problem of existence of NPPT-BE, which is still unanswered, is closely related to the 
Hubert's 17th problem. In Sec.|II] we cast the original problem into a 2 x 2 block matrix form. In Sec.|IIIl we find the relation 
between the original problem and the Hilbert's 17th problem. A similar but slightly different approach to the problem is also 
described. Section|lV]is conclusions. 

n. 2x2 BLOCK MATRIX FORMULATION 

To begin with we write a Schmidt rank-2 state vector in •7/®'^ i8> 'H^'^ as \i/f2 ) - Zi/t=i,2 Q |(Ai) with 

d d 

1*^*^ " Yj Yj 4l..i/hj2...h I'1'2 ■ ■ • '^>^ ® l^'l^2 ■ • ■ JN^B . 

i]j2,--,iN = l JiJi,--Jn=1 

where y!j._i tC? = 1, Yf , ,■ , u . u . - Sk,h, and Y/; , , , v. ' . v . = 6k,k^ due to the orthonormal condition. 

Our goal is to show {<{^P\ K{a)^'^ [/z^^) > for l/d < a < 1/2 and for all [i^f ^), where K{a) = / - daP+ is defined in Eq. ©. 
We write 

d d I V (aA^^ mI'V] \ / \ 

with Wi; = (qw, / , , CkU. , o, . . . , CbU.. .Y . In Eq. Oil, matrix elements of M), , , are bilinear forms of v. . . and v. . and 

■""(l,!) (1,2)- 



For A^ = 1, it is easy to see M\.' , = Vl"A{a)Vl', where 



(ki,h) h 



ym^ 



(y(k) 

y^'^ 
... 



!,(« 



(4) 



isad-^ xd matrix with v^*^^ = (v\ ,V2 ,■ ■ ■ ,vj)^,hy which one finds 

V (k,,k2)Jiii2...iN.j]J2--JN 

d d 

This is the coefficient of u''^'^* u'-'"^* . . . m*'^*m<'=^^m*'^^' . . . m**^^^ in (Sk, I K{aY'^ W) with 



(5) 



I0*> = 



\i,M 



Hence, if we replace v. v. ... v to v. ? . in vl and read u u . . .u us u . , Eq. (|5]l gives the coefficient of 

' -■- 'l '2 ''>' l\l->...lM k h l-> Im 1\Ii...Im^ ■* ^——* '^ 



'N I\l2---'N 



u- j^* ■ u. -^ . in {ij/ki I A(Qr)®^ I'Ajti)- Consequently, we can write 



'1 '2 'A' '1'2---'A' 



where y^ ' is a li^'^ x c/^ matrix obtained from y| '® by the above prescription. By this construction of the Hnear map, 



y[Ar] ^ ■^ ^(i)( .^ ^ .^^ _ ^ ^ .^^^. ^ y .^ ^ . . . ^ y.^ 



(6) 



!l,!2, ...,/«=! 



with 2f ; ; ^1 |f**'(;i,;2, ■ ■ -^'Af)! = 1, according to the expansion of the normaUzed vector v**^^: 

d 

y W ^ ^ ^W(;j ^ /2, . . . , /a-)v„ ® V,, ® ■ ■ ■ ® V,„ , 

!l,!2, ...,/>=! 

where {vi, V2, . . . , v^) is a set of orthonormal basis of the J-dimensional complex vector space C''. In Eq. dU, y- is given by 
Eq. dUi with V**' -^ v,-. 

It is convenient to expand a (normalized) vector \x\. ) in C' on which the operator mL ^ , acts as follows. 



kr)= Yj rf''^(iuh,---,iN)\xi,)®\xi,}®---®\xi^). 



i\ ,/2,.--)'A' — 1 



Here, Yfi i ,- =i \^''''\h,h, ■ ■ ■ , in)] - 1 and |jc,) = v* (/ = 1, 2, . . . , li), which also constitute a set of orthonormal basis of C''. 
This choice of |x,) is due to the observation that v* is the eigenvector of y. A(a)y, with the eigenvalue I - a. This fact is easily 
verified by direct computations. Writing the vector Vj |x,) explicitly as 



^.■v; = 


{\j . 
V,. 


■ °1 


((viTA 


_ 






.0 ... . 


.. vJ 


■(^'O:^ 




l(v,)*v,J 



the inner product between this and the vector \pp) is readily computed as {pp\ V^v* - (yiy{vj)p for p - 1,2,. . .,d,hy which 
one can show that 

d d 

{xi\ v]lV{a)Vi\x,„) = (v;)"'"y]'A(Q')y/v,; = vfyjV/v,; - a^(v,)p(vj)* ^(v/)^(v„,)* = 61^6 ji - a6ij6i,„. 

p=l q=\ 



In the following, we investigate the positivity of (i/'j | A(a)®'^ P2 )' ^■^■ 



■'"(a,!) ■'"(2,2)/ 



> 0, (7) 



or equivalently 



fL«Jr„-v*" v^^> w^" w<^>V- (1.1) ™(1,2^ 

/^ (Qf,v ,v ,w ,w ).-| (2) \j^m j^m 

I V™(2,I) ■'"(2,2)^ 



w(i) 



>0 (8) 



for all v('>, v<2', w('\ and w^^^ when 1/^ < a < 1/2. 



m. SEARCH FOR POSITIVE SEMI-DEFINITENESS OF M'"! (a) AND THE HILBERT'S I7TH PROBLEM 

If the polynomial /j on v*'\ v'^^\ w*'\ and w'^^ defined by Eq. (O were written as the sum of squares of the absolute values 
of polynomials, it would be positive semi-definite. However, this is not the case even for A^ = 1 in which /j "Ms known to be 
positive semi-definite when 1/d < a < 1/2 ll20ll22r . To see this we assume that all variables involved are real and set li = 3. 
Furthermore, we set jCj - x\^ - 3'3 - 3'3 - z&'R for simplicity. The polynomial fj^^^ \a - I /2) is now computed as 

/]!3='](a = l/2;v(i>,v(2),w('\w(2)) =/(v<'\v('>,vf vf ,w<'\vv<'\wf vvf ,z) 

. 2,^ + ^2^(;)2 + ,2,(')2 + 2z2v(;>vf> + z2,(2)2 ^ 2,2,(1),(2) + ^2^,(2)2 _ 2,2^(0^(1) + ^2^(1)2 ^ ^,(1)2^(1)2 ^ ^(1)2^(1)2 

T 2 (1) (1) (1) (1) (1) (1) , 2 (1)2 , (1)2 (1)2 , 1 (1)2 (1)2 ^ 2 (2) (2) , ^ 2 (1) (2) , (1) (2) (1) (2) 

, T (1) (2) (1) (2) (1) (2) (1) (2) , 2 (2)2 , 1 (2)2 (2)2 , (2)2 (2)2 ^ 2 (2) (2) (1) (2) (1) (2) , ^ 2 (1) (2) 

H- 2v\^ Vj '^i '^i ~ ^2 M ^^2 "^1 "*" -^ '^l "*" o M "^i + ^2 '^l ~ 2-^ "^2 "^2 ~ M ^2 "^1 ^^2 "*" ^-^ ^^2 ^^2 

, T (1) (2) (I) (2) , (I) (2) (1) (2) (2) (2) (2) (2) , 2 (2)2 , (2)2 (2)2 , 1 (2)2 (2)2 .^.s 

H- 2Vj Vj Wj Wj H- Vj Vj Wj Wj " Vj Vj Wj Wj "H Z Wj "H V j Wj "H - Vj Wj (9) 

This is the homogeneous real polynomial of n^ar - 9 real variables with even degree 2df - 4. 

We follow the relaxation method |43] to denote the column vector whose entries are all the monomials in v* , Vj , . . . , Wj , 
and z of degree at most df - 2 as X. Furthermore, let X./ denote the set of all real symmetric « x « matrices M(a = 1 /2) such 
that / - X^M{a - 1 /2)X, where n is given by 

If there exists a positive semi-definite matrix M{a - 1 /2) > within ths set X/, / is written as the sum of squares of polynomials 
(SOS form). Our task is to find such a positive semi-definite matrix. Note that the monomial entries of X with degree other than 
two can be omitted from the start because the squares of such components do not contribute to the homogeneous polynomial / 
of degree four; the corresponding diagonal elements of M{a - 1 /2) are all zero. Thus, the vector X is reduced to the following. 

V f 2 (1) (1) (2) (2) (1) (1) (2) (2) 

X - l^z , zv\ ',zv\ ,zv\ ', zv\ ',zw\ ',zw\ ',zw\ ,zw\ , 

(1) (1) (1) (1) (1) (1) (1) (1) (2) (2) (2) (2) (2) (2) (2) (2)1 

Vj Wj , Vj Wj , Vj Wj , Vj Wj , Vj Wj , Vj Wj , Vj Wj , Vj Wj | , 

The reduced set Hf consists of all 17 x 17 matrices with c,j- e R as follows. 

1 \ 1 {Ma Mc\ 



'^^"^2r2\Ml M, 



where 



with 



Ma = 



fi,i 






M'/- 



c4 o^ 

2 2 

2 2 

2 2 

2 2, 



Mr- 



/'2 2 0^ 

2 2 

2 2 

2 2j 



and 



< = 





-2 - cijo 

-Cl,l2 








-ci,ii 

-2-C1J3 











-2-ci,i4 

-Cl,l6 







-Cl,l5 

-2-ci,i7y 



with 



M„ = 



fi,i 






/' 






1 



ClO,13 





2 

-1 -Cio,l3 










cm,n 


- c 


0,13 





2 












1 



IVlg 



1 




Cl4,17 





2 

-1 -Ci4,17 










Cl4,17' 


- c 


14,17 





2 












1 



and 



M 



1,2 



(I -n 

2 
2 
-10 1 



and 



Mr 



ci,io 


Cl,ll 


Cl,12 


Cl,13 


Cl,14 


Cl,15 


Cl,16 


^1,17^ 








C2,12 


C2,13 














-C2,12 


-C2,13 






































C4,16 


C4,17 














-C4,16 


-C4,17 











C6,ll 





C6,13 














-C6,ll 





-C6,13 
































C8,15 





C8,17 














-C8,15 





-C8,17 


j 



One can check the identity / = X^ M{a - 1 /2)X by direct computations. To investigate the positivity of M{a - 1 /2), we make 
following replacements, ci lo — > ci, ciji — > C2, . . . , c^ n — > c\^, for simplicity. Here we denote the principal submatrix (PSM) 
that lies in the rows and columns of M(q' - 1/2) indexed by i,j,k as M^ij^^. We recall here that the necessary and sufficient 



TABLE I: Values of c, 



PSM 


variable 


PSM 


variable 


■'1^12,6,81 


C\ 


= -2 


^{3.11.15) 


cio = 


^(2,7.91 


C2 


= 


A/{4,12,16) 


cii=0 


^3,6.81 


C3 


= 


A?|5, 11,151 


Cl2=0 


^(3,7.91 


C4 


= -2 


'^16.11.151 


Cl3=0 


^{4,6.81 


C5 


= -2 


'^17,12,161 


Cl4 = 


■'1^(4,7,91 


Cf, 


= 


^18,1 1,15) 


Cl5 = 


^(5,6.81 


c? 


= 


A?{i),12,16( 


Cl6 = 


^(5,7.91 


Ci 


= -2 


A?|l 1.12.161 


Cl7 = -1 


W|1.12.16| 


C, 


= 


M|i2,15,16| 


Cl8 = -1 



condition of M{a - 1/2) > is that all the PSM's of M{a - 111) are positive semi-definite |142|]. The PSM M(2,6,8) takes the 
form 



( 



M, 



|2,6,8) - 



-2-Ci 





-2-ci 0^ 
2 2 
2 2, 



Requiring that all the eigenvalues of M|2,6,8| be non-negative, we have ci = -2. Similar procedures determine remaining variables 
from Q to ci8. The results are listed in Table Hill 

Using these values of Cj (i - 1,2,..., 18), the minimal eigebvalue of M(a - 1/2) is found to be 1 - Vs < 0, which means 
that positive semi-definite polynomial / does not take the SOS form. 

For a = 1 /<i = 1 /3, we also define the matrix M{a - 1 /3) such that f^'lj^Ha - 1 /3) = X'' M{a - 1 /3)X. Similar procedures 
described above lead to the following form for M(a - 1/3); 

^ 1\ ifM, Mc^ 



3 / 3 {Ml Mb) 



where 



Ma 



f8 0000000 0'' 
030300000 
003030000 
030300000 
003030000 
000003030 
000000303 
000003030 
0000030 3, 



M« = 



(2 0-1 
3 
3 
-10 2 
2 0-1 
3 
3 
-10 2 
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2 j 



and 



Mr 
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0. 



The minimum eigenvalue of M{a - 1/3) is found to be zero. Since /j " '(a) is linear with respect to a, the maximum of 
^mm{M{a)), the minimum eigenvalue of M{a), with respect to ci {i - 1, 2, . . . , 18) is negative for l/d < a < 1/2, the positive 
semi-definite polynomial /j^^ '(a) cannot be expressed as an SOS form for 1/3 < a < 1/2. 

The problem of existence of such positive semi-definite polynomials that cannot be written as an SOS form is related to the 
Hubert's 17th problem that dates back to 1900. The Hilbert's 17th problem was solved affirmatively by Artin in 1927. 

Theorem (Artin) — Every polynomial p e R[x] that is nonnegative on R" is a sum of squares of rational functions, i.e., 
P - Hji^j/bj)^ for some Uj, bj g R[x]. 

Here, one may ask a question. Is it possible to write every positive semi-definite polynomial as a sum of squares of polyno- 
mials? Unfortunately this is not true. A series of counterexamples to this question was firstly found in 1965 by Motzkin. One of 
the special cases of Motzkin's polynomials takes the form 

Pm - XjX2(xj + X2-2>) + I, 



which is nonnegative but cannot be expressed as an SOS form. Since then many such polynomials are found. Our polynomial 
/ [Eq. ©I consists one of them. See |44] for an excellent survey on the Hilbert's 17th problem and the historical remarks 
including Motzkin's polynomials. To show the positivity of/, we have to construct a set of rational functions satisfying Artin's 
theorem. This is actually quite a formidable task. However, Reznik proved the following recently [45]. 

Theorem (Reznik) — Let p e R[x] be a positive homogeneous polynomial on R". \f p > on R \ {0), then there exists r e N 
for which the polynomial (2;Li x^) /? is a sum of squares. 

This theorem might be helpful to prove /j > 0. However, small r trial does not seem to work well even for A^ = 1 and real 
variables case. What is the minimum r - r(N, d) to show /] ' > 0? Is it impossible for general A^? These questions remains to 
be answered. 

So far we have utilized Eq. ^ to investigate the positivity of M^^\a). This is not an only option. The inequality [Eq. (I8]l] 

is equivalent to the following M- (A) M^^\^ > and (B) (^^^1 ^lu) VT) (4'^1 ^1^2) VT) ^ K-^T 



^(U)kr)| for all 



The first inequality is proven as follows. We decompose A(a) into Ma) = (1 - a)I + aZ so that 



ix; ^) 



.fJ^ 



(x['^i|y[^i'A(a)«^y['^]|x['^i) 

N 

= ^(i-<a^-"((x[^i|y|'^"7« 



■»{N-n)^AN] \\N] 



>Y. 



M 



>..^^(,, 



[N]\ ^AN\] y<s{N-n) , 



v 



vT' wr')) ■ 



.imym^^m 



(10) 



Here the special term {x\ 



[N]\,AN]i^m,AN] \JN] 



1 



y|"J'z«"y;-'U' 



> is given by the following sum of squares. 



'\<]\ 



lN<jN 



" Z • ■ ■ Z Z ■ • ■ Z 'S"^'^'^ • ■ ■ ssn(PN)^'\Pi(iuji), . . . , PnUn, JN))r]^'^*(PiUuii), ■■■, PnUn, iN)) 



(11) 



where Pk{aii,bk) is a^ (sgn(Pi) = 1) or b^ (sgn(Pk) = -1). Terms other than this are obtained by replacing some (or one) Z's 
in Z®^ by /. For example, when the p-th Z is replaced by /, the corresponding SOS form is obtained by dropping sgn{Pp), 
replacing Y.i,<j„ by Z;, , Ijj,,^ PpUp, jp) in ^^'^ by i,,, and Ppiip, j,,) in rf^'^ by j,, in the right-hand side of Eq. (??). One can prove 
the above results by induction. Since all terms in the right-hand side of Eq. (fTOl l are (strictly) positive for l/d < a < 1 /2, we 
conclude that m\^^ > 0. 



As for the second positivity condition (B), it is very hard to find a closed form for 

^ Vl ^"(l.Dri /Vl I™(2,2)r2 / |\-^l I ■'"(1, 2) 1-^2 /| ■ 

However, if we assume that all variables are real, we can prove the following by induction, 

d 



Q[N=i] = i J] J] J ([ijkk] - [jikk] + [kkij] - [kkji] + [kijk] - [ikkj]) 
(=1 j=l [k=l 

-. d d d d 
^JrYjYjYjYj [^^1 ~ ^3 + gsf + igl -gA+ gbf + {gl -gi+ gbf + igl -gA+ gif\ : 



!=I j=\ k=\ /=1 

where [ijkl] is the abbreviation of v- v^- w\ w\ , 

g\ = [ijkl] - [ijlk] - [jikl] + [jilk], 

g2 = [klij] - [klji] - [Ikij] + [Ikji], 

g3 = [ikjl] - [iklj] - [kijl] + [kilj], 

g4 = [jlik] - [jlki] - [Ijik] + [Ijki], 

g5 = [iljk] - [ilkj] - [lijk] + [likj], 
and 

g6 = [jkil] - [jkli] - [kjil] + [kjli]. 

It seems to be very hard to generalize this result to N > 2 and complex variables case. The possibility of non SOS character or 
nonpositivity of &^^ cannot be excluded. 

IV. CONCLUSIONS 

A trial to solve the problem of existence of NPPT-BE has revealed the intriguing relation to the Hilbert's 17th problem. We 
have witnessed much progress of commutative algebra and real algebraic geometry concerning the Hilbert's 17th problem such 
as Reznik's theorem recently. One may expect that some of such concepts and techniques developed so far eventually lead to the 
final proof or disproof of the existence of NPPT-BE. 
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